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Abstract 

Jahn- Teller (JT) electron-phonon coupling effects in the colossal magnetore- 
sistance perovskite compounds La\- x A x MnOj, are investigated. Electron- 
electron correlations between two degenerate Mn e g orbitals are studied in 
the Gutzwiller approximation. The static JT distortion and antiadiabatic po- 
laron effects are studied in a modified Lang-Firsov approximation. We find 
that (i) the electron or hole character of the charge carrier depends on the 
static JT distortion, and (ii) due to the two-component nature of the JT 
coupling, fluctuations in the JT distortion direction contribute to the charge 

transport in similar fashion as the local spins. 
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The recently rediscovered "colossal magnetoresistance" (CMR) perovskite-based 
Lai_ x A x MnOs compounds (A = Ca, Sr, Ba) exhibit striking magnetic field and temper- 
ature dependence of the resist ivityS^. Strong correlation between the electrical resistivity 
and ferromagnetism in these materialsi was qualitatively explained by Zener's double ex- 
change (DE) modeli. In this model, the three t2 g Mn d xy , d xz , d yz orbitals are completely 
filled and the electrons are localized; thus their spins can be treated as a localized com- 
posite spin S = 3/2; the conduction electrons are provided by e g electrons on the Mn 3+ 
ions. Strong Hund's-rule coupling of strength J H 3> t (the hopping matrix element) of 
the electronically active e g electrons to the localized t2 g electrons, results in the motion of 
electrons spin-polarizing the localized Mn spins and leading to ferromagnetism. An applied 
magnetic field will polarize the local spins and thus decrease the electrical resistivity. The 
DE model was subsequently refined by Anderson and HasegawaB, and de Gennesi, and 
studied more recenthJHll to explain topical experiments. However, it now appears that the 
DE model can not adequately explain various experimental data by itself: Jahn- Teller (JT) 
electron-phonon coupling also plays an essential 

Previously, we studied a model of combined DE and JT coupling^ and found that the 
JT coupling can substantially reduce the resistivity and the Curie temperature T c . Using 
this model, we can qualitatively understand the metal-insulator (MI) transition around T c : 
namely, a polaron localized by the combination of the changing of the spin polarization cloud 
around it and Anderson localization due to the dramatic increase of spin-disorder scattering 
near T c . Here we study the DE model with JT coupling using two e g orbitals at each lattice 
site. One motivation is to incorporate the doubly-degenerate states to take into account 
the mobile carrier density correctly and study the electron-electron (e-e) on-site correlation 
effects. A second motivation is that the doubly-degenerate states act like a pseudospin which 
couples to a vector (of two phonon modes). For T close to or above T c , a pseudospin-induced 
"double exchange" effect will then become important (see below). 

The Hamiltonian we use consists of the DE part 7i DE and coupling to the lattice repre- 
sented by H JT . H DE can be derived from the Kondo-lattice Hamiltonian 



(i,j)a i 

— Jh^Z Si ■ <7 a p(a ia aip + b ia bi^i), (1) 

i 

where a is the Pauli matrix and Si is the local spin of the three t 2g electrons (S = 3/2). The 
operators a ia (a\ a ) and b ia (b\ a ) annihilate (create) an electron at the two e g orbitals (a and 
b) with spin a at site i. The second term in Eq.(|l]) describes the e-e interaction between 
a and b orbitals on the same site. A detailed calculationH showed that for the Hund's-rule 
coupling Jh 3> to/2, the eigenstates of a two-site systems are approximately formed by 
electrons all parallel (or all anti-parallel) to the on-site local spins. Based on the counting 
of single electronic states, it is natural to postulate this will also be true for iV-site systems. 
So in the limit of Jh 3> to/2, Ji DE can be written as 

H DE = -to E Tlfe + btbj + h.c.) + Vj2n?nl (2) 

(iJ) i 

where we have omitted the spin index of the electrons since the electronic spin is fixed by the 
local spins: the single-electron states form two bands with a gap ~ Uu ^ ksT. Thus only 
the lower band is occupied, and it can be assumed there is no double degeneracy of each e g 
orbital. (This assumption is further strengthened if there is a large on-site Hubbard inter- 
action). 7s is the bandwidth renormalization due to the DE mechanisml: 7^ = / s ° 2 g^ \ , 
where Sq is the total spin of the subsystem of the two localized spins on sites i and j and 
the electronic spin. The angular brackets < . . . > denote a thermal average. 
H JT can be written as 

i 

+zf(^+3y+£f (qi+qi), (3) 

i i 

where Qi = (l/y/6)(26Z - SX - 5Y) and Q 2 = (1/V2)(5X - 5Y) are the two local JT 
distortions0'0, with SX 1 the distance change between neighboring 2 ~ ions along the X 1 - 
axis. Here we have used pseudospin notation to represent the two degenerate orbitals: 
= a; cj = b with a denoting the pseudospin index. Then r| = o^q^qs. When we define the 



— * — * 

distortion vector Q by Q x = Qi and Q z = Q 2 , the JT coupling becomes XjtQ-t. This shows 
clearly that the eigenenergy of the system is independent of the distortion direction Q/\Q\. 
If we assume that \Q\ is fixed, Q can be treated as a local "classical" spin. Then the JT 
coupling is formally similar to the electron-spin coupling in Eq. ([]]). In the CMR perovskites, 
if the JT electron-lattice coupling is close to the adiabatic limit and the coupling is strong (i.e. 
effective \jt\Q\ ^ to/2), it is thus possible effectively to have additional "double exchange" 
from the 2-component-doublet JT coupling. However, there are several mechanisms which 
will stablize a specific distortion direction, including (1) anharmonic lattice interaction^, 
(2) higher-order JT interactions^, (3) interactions between local distortions at different 
sites0. To be specific, we consider only effect (3) here, by introducing an effective near- 
neighbor "ferromagnetic" coupling. At low T, the JT distortion will be homogeneous in 
a certain direction. However, at higher T, fluctuations of the local distortions will induce 
the additional "double exchange" effects. The local JT distortion will start to be observed 
for T < rjl. Since Q is a classical variable constrained to the x — z plane instead of a 
sphere, the fluctuations of Q/\Q\ will be larger than a corresponding local spin S. Thus 
its fluctuation is still significant at T < T c even though the Curie temperature T c is smaller 
than the "paramagnetic transition temperature" of Q/\Q\ alone. 

The following calculations are restricted to the metallic phase (T < T c ). Incorporating 
the above JT distortion-induced "DE" effect introduces an additional bandwidth reduction 
factor 7jt. At low T, the distortion direction is "frozen", 'Jjt = 1. At T closer to T c , 
we will approximate the stablizing force by a "ferromagnetic" coupling between the Q^s. 
The strength of this coupling is unknown and needs to be determined by (e.g. structural) 
experiments. After a rotation to make the distortion direction point along the z-axis and 
quantization of 7i JJ ', the Hamiltonian components become 

H JT = -fipto E maO&Bl + B t ) + huJ2 BjBi, 

ia i 

n DE = -t E ih%(4^ + h.c.) + vJ2 n im . (4) 

(i,j)a i 

Here, the electron-phonon coupling strength is e p = (Xjt) 2 /2K , and fku is the phonon 
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frequency. 

We treat the dynamical phonons within the homogeneously modified variational Lang- 
Firsov approximation^ using a canonical transformation IA = e~ Sl ^ A ^e~ S2 ^\ where 
Si = 1/2 J € p huj J2i A(B}—Bi) is designed to describe static lattice distortions through the in- 



troduction of a (variational) static displacement field A. £2 = —\J ^ P /^Y^ia1 n iaC^(Bl — Bi) 
describes the antiadiabatic polaron formation, with 7 measuring the degree of the po- 
laron effect. After the transformation U we can average 7i over a squeezed phonon state 
\^ ph ) = exp(-EjT(S i S i - B}Bl))\vac) with r a variational paramete r§. The e-e corre- 
lations are approximated by the Gutzwiller approximation^. The free energy can then be 
written as 

T = -X a (Pa + d)- X b ( Pb + d) + (V + 2e p7 (2 - 7 ))d 
~\ E ln[l + e-^^-^)] -e p (2 7 - 7 2 )(l-x) 

P k,i£{a,b] 

+ + r- 2 ) + + (1 - 7 )A(p a - p b ) +F s + F q (5) 

with the electron dispersion relation given by = 2 7s 7jt^o^(t) (cos^) + cos{k y ) + cos(A^)) 
with a polaronic band narrowing factor £(7) = exp(— e P 7 2 r 2 /(/iu;)). Here we have used 
the notation f — > a; | — > b. are the discontinuities in the single-particle occupation 
number at the Fermi surfaces, which give additional band narrowing due to e-e interactions: 
m*/m ~ 1/qi- The relations between and Pi (probability of occupation of orbital i) and 
d (double occupancy) is 

qi ~ (p i + d)(l- Pi -d) ' 

where a = b and b = a. x is the doping density. Note the crucial distinction between 
the static JT and antiadiabatic polaron distortions: the double degeneracy of e g is lifted 
by (1 — 7) A; in contrast, the antiadiabatic polaron distortion contributes an effective e-e 
interaction 27(2 — 7)e p . 

T s and T q in Eq.(|5]) are free energies from the local spin 5 = 3/2 and JT distortion 
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direction Q/\Q\, respectively. We used molecular crystal approximations for both JfJH and 
T q . T q contributes only at high T: 



T q = leff{Q Z ) ~ \0g{v JT )/(3 - 6J ef f(Qz) 2 - 



(7) 



Here l e ff is the effective field for Q/\Q\, and v JT = I (j3l e ff). (Q z ) = Ii((3l e f f) / I$(f3l e f f) . In 
this molecular crystal approximation, ^ JT = sinh(2/5/ e //)/ (7r/3i e ///Q (Pl e ff)), where I n (x) is 
the n-th order modified Bessel function. 

The mean-field equations at low T {^jt = 1) are derived by minimizing the free-energy 
Eq.(§: 



ie{a,b} J °Pi 



1 — X 



J2 J dep{e)f{q h \ h e 



i€{a,b} 



i€{a,b} 



(V + 2e p7 (2- 7 ) 
(l- 1 )[(l-x)-(p a - Pb ) 2 



r 2 r 

7^ X, y rfe p( e )(^ e )/(^> 



46 p7 2 
(M 2 



ie{a,6} 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



Here / A», e) = [1 + e /3 ^ e_/i+Aj ' ) ] _1 is the Fermi-distribution function, and p(e) is the 
electronic density of states. For the dispersion relation = 2t e f f(cos(k x ) + cos(k y ) + cos(k z )) , 
pie) = ^TJJ io° ^ x cos (sf 2- ) [^(^c)] 3 ! with Jo(^) the zeroth order Bessel function. Below, 
we evaluate = Jf^ dxp(x), which converges easily and can be accurately approximated 
with high order Chebyshev polynomials. 
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Eq.(|8|)-(|13|) are solved numerically by iteration. The resistivity is calculated using the 
same approximation as in Refl, which is comparable to the memory function method0. 

The resistivity is given by: p = p a Pb/(p a + Pb), where pi oc ^-^A(7 S 7 JT ) with m* ~ l/t eff 

f 1 

and 

A(7^) = ((|^) 2 )(cos 2 (%)) - { lslJT f (14) 

with 7j T = 1 at low T: A (7,7^) = ((^f^) 2 ) - (7 S ) 2 . At higher T, the fluctuation of 
local JT distortions begins contributing to the resistivity. The self-consistent equation for 
l e ff can be derived easily from minimizing the free energy, and is similar to Eq . (|13D through 
the substitution: h — > l e ff, m — > (cos |), and 7 S — > jjt- The resistivity is readily calculated 
given (cos 2 (#)> = 1(1 + If((3l eff )/P ((3l eff )). 

In Fig. [l|, we show the variations of the electron density of each e g channel. It can be seen 
that the character of the carrier is strikingly dependent on the static JT distortion A. For 
small A (or A = when there is no JT coupling), the charge carriers will be electrons, since 
the electron density in both channels is lower than half-filling. Only with large enough A does 
the carrier (in the majority channel) become hole-like. The e-e interaction V can change only 
the double occupancy but not the carrier character. (Nagaoka "ferromagnetism" is unlikely 
here due to the doping x and moderate V interaction strength). With large JT distortion, 
because of the creation of the majority channel, the e-e interaction V is unimportant for 
charge transport. Based on this result, if experiments (e.g. Hall transport or thermoelectric 
power) can detect the hole-character of the charge carriers, this will indirectly indicate the 
presence of a large static JT distortion, either local or homogeneous. It is indeed observed 
that the thermoelectric power is positive (hole-like) at x = 0.2 and negative (electron- 
like) at x — 0.4@, which is consistent with Fig. [I] in a certain range of electron-phonon 
coupling and phonon frequency. In our choice of model parameters we have used the hopping 
matrix element as the energy unit: the temperature unit is calculated with to = lei> , so if 
to = 0.2ev , then T = 400 will be T = 80K. The unit of magnetic field is renormalized by 
H = gp^sB, where g is the effective Lande (/-factor. Our self-consistent mean- field solution 
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gives T c ~ 2550 for x = 0.2, i.e. T c ~ 500K with t = 0.2ev. 

Fig. |2] shows the magnetic-field dependence of the magnetoresistance ratio Ap/p in 
the low-T (7jt = 1) region. The global features of the Ap/p — H curves are somewhat 
similar to the experimental data@. Note that we have not included any additional "residual" 
resistivity mechanisms (disorder, phonon scattering, etc) which will strongly decrease Ap/p 
at low T, and would need to be subtracted from experimental data for comparison with our 
prediction. From various experiments, Ap/po can be as large as 0.8 for moderate to high 
T around T c , so the spin-fluctuation scattering considered here will then dominate. The 
temperature dependence of the magnetoresistivity ratio Ap/po is shown in Fig. |3[ Without 
the local JT distortion direction fluctuation (j JT = 1), the Ap/p — T curve is flat close to 
T c . At high T (T <T C and jjt ^ 1), the effects due to the local JT distortion fluctuations 
can make | Ap/ po| increase with T, see Fig. |3], which has also been observecS Even without 
self-consistently including higher order JT coupling, our calculations qualitatively agree with 
the experiments: there is a characteristic temperature T p (T p ~ 1000 for the parameters we 
chose). For T < T p , the local JT distortion direction is frozen, the resistivity is then due to 
fluctuation of local spins and other "residual" scattering mechanisms. For T p < T < T c , the 
JT distortion direction fluctuation begins to contribute to the resistivity, and gives correct 
magnitude and T dependence of Ap/p . 

In conclusion, we have studied the combined effects of JT electron-lattice coupling and 
the double degeneracy of e g orbitals in the CMR perovskites in their metallic phase (T < T c ) . 
We found that without the static JT distortion, the carriers are electrons. However, with 
sufficiently large static JT distortion, the carrier can be hole-like. Also, the additional 
scattering effects of the JT distortion ( effective "double exchange" effects of the local JT 
distortion direction) will be enhanced in the inhomogeneous insulating phase (T > T c ) where 
small magnetic polarons are localized. This will be reported elsewhere. 
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FIGURES 



FIG. 1. Occupation density variation of each e g channel with the electron-lattice coupling e p 
at low T. to = 1, V = 6.0, hu = 0.5, T = 400. A increases monotonically with e p . 

FIG. 2. The variation of Ap/po with magnetic field at low T {^jt = 1). to = F e p = 2.0, 
F = 6.0, hu> = 0.5. 

FIG. 3. The variation of Ap/po with T at low (7jt = 1) and high T (T < T c ) with J e ff = 0.1. 
Parameters are the same as in Fig.2 and Ap = p(H = 0.03) — po- Note that the magnitude of 
Ap/pQ at low T can not be compared directly with experiments due to "residual" scattering (see 
text). 
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